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SECTION  1 


INTRODUCTION  TO  EQUATIONS  OF  MOTION  FOR  CONTROLLED 
SPINNING  PROJECTILES 

The  notations  and  conventions  which  have  been  adopted  are 
consistent  with  old  NACA  practice,  with  the  report  by  Charters, 
reference  [1]  ,  and  with  the  monograph  by  Jones,  referen.ee  [2]. 

Newton's  laws  for  describing  the  motion  of  rigid  bodies  state 
that  the  sum  of  all  external  forces  acting  on  a  body  equals  the  time 
rate  of  change  of  momentum,  and  the  sum  of  the  moments  of  forces 
equals  the  time  rate  of  change  of  moment  of  momentum. 


All  quantities  in  equations  (1-1)  are  specified  relative  to  axes  fixed 
in  space.  Consequently,  the  six  component  equations  of  motion  repre¬ 
sented  by  the  vector  expressions  of  equations  (1-1)  are  too  unwieldy 
for  practical  use,  since  the  moments  of  forces  and  the  body  inertia 
matrix  would  vary  with  time  with  respect  to  an  inertial  (fixed)  axis 
system.  The  difficulty  is  overcome  by  referring  all  quantities  to  an 
Eulerian  or  moving  axis  system  which  is  coincident  with  a  set  of 
axes  fixed  on  the  body.  The  formulation  in  terms  of  moving  axes  is 
made  by  recalling  that: 

•  -T-—  +  Ex 

fixed  axes  moving  axes 

» 

which  permits  equations  (1-1)  to  be  written  as 
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■'N 

zT  -  ■•d^jffi-  +  nx(mU) 

> 

zm  -  -dffi  +  fix(S) 

J 


(1-2) 


Equations  (1-2)  describe  body  motion  in  terms  of  linear  momentum 
and  moments  of  momentum  referred  to  the  moving  coordinate,  and  n, 
the  angular  velocity  vector  of  the  moving  frame  relative  to  inertial 
(fixed)  space. 


The  force  components  and  the  moment  components  are  resolved  onto 
a  coordinate  frame  which  is  associated  with  the  projectile. 

There  are  a  number  of  coordinate  frames  which  are  useful  in  the 
analysis  of  controlled  spinning  projectiles.  These  will  be  intro¬ 
duced  at  this  time,  even  though  most  of  them  will  not  be  used  until 
much  later.  The  Justification  for  introducing  them  here  is  that 
most  of  this  notational  complexity  is  caused  by  the  multiplicity  of 
coordinate  frames. 


Each  coordinate  frame  is  a  right-handed  coordinate  frame  with 
origin  0,  and  three  mutually  orthogonal  axes  which  are  designated 
the  x-axis,  y-axis,  and  z-axis,  respectively.  A  right-handed 
coordinate  frame  is  one  in  which  the  conventions  for  positive 
rotations  have  been  coordinated  with  the  positive  direction  of  the 
axes  so  that  a  right-handed  screw,  pointed  in  the  positive  direction 
of  an  axis  will  advance  in  a  positive  direction  when  rotated  in  a 
clockwise  sense,  as  seen  looking  in  the  positive  direction.  A  right- 
handed  frame  is  shown  in  Figure  (1-1). 
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Figure  (1-1).  A  Right-Handed  Coordinate  Frame 


The  following  coordinate  frames  will  be  useful. 

x°  y°  z°  an  inertial  coordinate  frame 
F  F  F 

x  y  z  a  coordinate  frame  fixed  to  the  cruciform  canards 

x'  y'  x'  a  coordinate  frame  fixed  in  the  body  of  the  projectil 

x"  y"  z"  a  stabilized  coordinate  frame  which  moves  with  the 
projectile 

T  T  T 

x  y  z  a  coordinate  frame  which  is  fixed  in  the  target 

x  y  z  an  aeroballistic  coordinate  frame 


Several  vectors  will  be  resolved  onto  the  several  coordinate 


frames  as  follows; 


R 

at 

ix 

+ 

jy 

+ 

kz  ; 

position 

U 

m 

iu 

+ 

jv 

+ 

kw  ; 

velocity 

A 

■ 

ia 

+ 

Jb 

+ 

kc  ; 

translational  acceleration 

ib 

■ 

ip 

+ 

jq 

+ 

kr  ; 

angular  velocity 

r" 

4b 

m 

iP 

+ 

Jq 

+ 

kr  ; 

angular  acceleration 

r 

m 

iX 

+ 

JY 

+ 

kz  ; 

force 

M 

m 

iL 

+ 

JM 

+ 

kN  ; 

moment 

where  i,  j,  and  k  are  unit  vectors  in  the  x,  y,  z  directions, 
respectively , 


(1-3) 


The  equations  of  motion  of  a  rigid  body  have  been  given  by 
Jones.  The  equations  are  expressed  with  respect  to  a  moving  coordinate 
frame  whose  origin  coincides  with  the  center  of  gravity  of  the  rigid 
body.  The  equations  will  be  written  initially  with  respect  to  an 
unprimed  coordinate  frame  and  the  semantics  of  the  equations  will 
tighten  up  in  the  discussion  which  ensues.  The  equations  are: 

m[u  -  vr  +  wq]  ■  X  +  XQ 

m[v  -  wp  +  ur]  ■  Y  + 

m[w  -  uq  +  vp]  «  Z  +  Zq 


)  Force  Equations  (1-4) 


h^  -  h2r  +  h3q  -  L 
h2  -  h3p  +  l)1r  -  M 
h3  -  h1q  +  h2p  ■  N 


)  Moment  Equations 


(1-5) 
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h1  *  Ap  -  Fq  -  Er 

hn  *  Bq  -  Dr  -  Fp 

feg  ■  Cr  -  Ep  -  Dq 

where 

XG>  Yq,  Zg  are  components  of  a  force  vector  due 
to  gravity 

^l’  h2  ’  h3  are  comP°nen,ts  an  angular  momentum  vector 

A,B,C,D,E,F  are  components  of  an  inertial  tensor,  and 

A,  B,  C  are  moments  of  inertia 

D,  E,  F  are  products  of  inertia 

The  application  of  equations  (1-4),  (1-5),  and  (1-6)  imposes 
certain  practical  constraints.  First,  it  is  desirable  to  choose  the 
coordinate  frame  so  that  the  components  of  the  inertial  tensor  (A,  B, 
C,  D,  E,  F)  are  constant.  This  can  be  accomplished  for  a  projectile 
of  constant  mass  by  causing  the  coordinate  frame  "to  be  fixed  in  the 
body  of  the  projectile.  This  is  the  primed  coordinate  frame  x',  y * , 
x').  Writing  equation  (1-4)  relative  to  the  grimed  coordinate  frame, 
gives 


m[u'  -  v '  r '  +  w'q']  -  X'  + 

m[v'  -  w '  p '  +  u'r']  -  Y'  +  Yg 

m  [w '  -  u  '  q  '  +  v ' p ' ]  -Z'  +  Z'Q 


Force  Equations 
Relative  to 
Body  Frame 


(1-7) 
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The  angular  momentum  equations,  (1-6),  become 


A '  -F '  -E ' 


-F1  B'  -D1 


d-8) 


|_4J  L‘E'  -D’  c,j  [r,j 

where  the  coefficients  A'.  B' ,  C' ,  D' ,  E’,  F'  are  constants. 
Differentiating  (1-8)  gives 

r  h j  i  -  r  a '  -f  -  -e •  i  r  p '  i 


B’  -D1 


(1-9) 


-E '  -D1 


Let  the  matrix  be  inverted  so  that 


-F'  -E'  ;  E 

B’  -D' 


A '  -F '  -E 1 

-F '  B '  -D ' 

_E '  -D '  C ' 


(1-10) 


Then , 


E'1  [&i 


(1-11) 
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Combining  (1-11)  and  (1-5) 


• 

p' 

-  E"1 

L’ 

+ 

h2r’  “  h3q’ 

4 

q' 

M' 

+ 

h^P ’  -  h^r' 

r' 

N’ 

+ 

hjq’  -  hgP ’ 

Introducing  the  notation 

a’  -  X’/tn  ' 
b’  -  Y'/m  • 
c '  ■  Z  * /m 

, 

The  equations  are: 


u* 

■  v' r '  -  w'q 

’  +  a'  +  (X^/m) 

V’ 

■  w'p'  -  u'r 

'  +  b'  +  (Y£/b) 

w’ 

■  u'q'  -  v'p 

'  +  c'  +  (Zq/bq) 

V 

-  E-1 

L'  +  h£r’  -  h^q ' 

9 

q' 

M'  +  hgP'  -  h^r' 

r  ’ 

N’  +  h^q '  -  h^p ' 

(1-12) 


(1-13) 


(1-14) 


(1-15) 


Equations  (1-13),  (1-14),  and  (1-15)  are  the  most  general  nonlinear 
equations  of  motion  which  will  be  considered. 

Second,  this  most  general  case  can  be  constrained  by  imposing  a 
rotational  symmetry  about  the  x'  axis.  The  full  significance  of  this 
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symmetry  can  best  be  implied  in  terms  of  the  inertial  coefficients. 
The  inertial  coefficients  are  defined  as  follows  (see  reference  [3], 


p.  139). 

A' 

3 

y'2  +  z'2)pdT 

*/T 

B' 

- 

/( z'2  +  X' 2)pdT 
*'t 

C' 

- 

Ax'2  +  y '  2 )  pdT 

D' 

- 

/ y'z'  pdT 

E’ 

- 

/ z'x'  PdT 
*> 

F’ 

- 

Jrx'  y'  PdT 

4 

where 

p  “  density 

dr  •  volume 

increment 

Lot  the  x',y',z'  axes  (the  primed  coordinate  frame)  be  aligned 
initially  to  coincide  with  the  x,y,z  axes.  Then,  let  the  x'  axis 
and  the  x  axis  remain  coincident  while  the  primed  coordinate  frame  is 
spun  about  the  coincident  axes.  The  inertial  components  about  the 
unprimed  axes  will  generally  not  be  constant.  However,  they  will 
be  constant  provided  the  coincident  spin  axes;  i.e.,x’  axle  and 
x  axis,  are  axes  of  summetry.  Under  these  conditions 


A  is  constant 
B  -  C 

D  *  E  ■  F  ■  0 


(1-17) 


When  these  axial  symmetry  conditions  are  satisfied,  the  equations  of 
motion  can  be  written  in  the  following  form: 


.  *N 

u  -  vr  -  wq  +  a  +  (Xr/m) 

4  W 

v  ■  wp  -  ur  +  b  +  (YQ/m)  \ 

w  -  uq  -  vp  +  c  +  (ZQ/m) 


(1-18) 


1-8 


(1-19) 


h1  -  Ap' 
h2  -  Bq 
h3  -  Br 

P'  ■  (L/A) 
q  -  (M  +  h3P 
r  -  (N  +  h1q 


J 


V)/®  > 

h2P)/B 


(1-20) 


Note  that  the  angular  momentum  component  (h^)  depends  on  the  spin  rate 
of  the  body  (p')  and  is  independent  of  the  spin  rate  of  the  unprimed 
coordinate  frame  (p). 


The  usefulness  of  equations  (1-18)  through  (1-20)  depends  on 
whether  the  force  components  (X,  Y,  Z)  and  the  moment  components 
(L,  M,  N)  can  be  expressed  conveniently  in  the  unprimed  coordinate 
frame.  They  can  be  and  the  manner  of  doing  so  has  been  described 
by  Charters  [1],  pp.  27-33.  According  to  Charters: 

"...  It  may  be  concluded,  therefore,  that  the  aero¬ 
dynamic  coefficients  are  invariant  with  respect  to 
rotation  of  coordinate  axes  about  the  axis  of 
symmetry.  Consequently,  the  y  and  z  axes  may  be 
oriented  at  will  around  the  x  axis,  without  regard 
to  the  orientation  of  the  missile  about  x,  since  the 
aerodynamic  coefficients  are  solely  functions  of  the 
missile's  external  shape  (and  such  non-dimensional  para¬ 
meters  as  Reynolds  and  Mach  numbers)  and  it  has  been 
shown  that  the  aerodynamic  coefficients  do  not  change 
with  orientation  of  the  missile  in  roll.  In  fact, 
the  missile  may  be  allowed  to  spin  about  the  x  axis 
with  respfect  to  the  y,a  axes  and  the  aerodynamic  co¬ 
efficients  will  be  unaffected." 
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Charters  defines  three  sets  of  coordinate  axes.  One  set,  which 
he  calls  "earth  axes"  correspond  to  the  inertial  coordinate  frame, 
xc ,  y°,  z° .  A  set,  which  he  calls  body  axes,  corresponds  to  our 
coordinate  frame  fixed  in  the  body,  x',  y',  z'.  A  set,  which  he  calls 
"pseudo-stability  axes"  corresponds  thus  far  to  the  as  yet  unnamed 
unprimed  coordinate  frame,  x,  y,  z.  Charters  points  out  that  the 
pseudo-stability  axes  may  be  oriented  at  will  about  the  axis  of 
rotational  symmetry.  The  pseudo-stability  axes  are  commonly  referred 
to  as  aeroballistic  axes. 

It  should  be  pointed  out  that  Charters  pseudo-stability  axes 
are  not  associated  with  a  physical  body.  Thus,  the  magnitude  of  the 
angular  velocity  about  the  x  axis  can  be  specified  arbitrarily  and 
does  not  have  to  satisfy  an  equation  of  dynamic  equilibrium. 

A  coordinate  frame  will  be  associated  with  the  canard  assembly, 
or  fin-frame,  of  a  controlled  spinning  projectile.  This  will 
require  an  equation  of  dynamic  equilibrium  which  is  introduced  in  the 
following  paragraphs. 

The  geometrical  relationship  between  the  primed  and  unprimed 
coordinate  frames  is  the  same  whether  or  not  the  unprimed  frame  is 
physical.  This  relationship  is  shown  in  Figure  (1-2). 

One  useful  treatment  of  the  unprimed  coordinate  frame  is  to 
set  its  spin  rate  (p)  equal  to  zero.  The  motion  of  the  projectile 
is  then  described  relative  to  a  non-spinning  coordinate  frame  and  the 
equations  are 


u  ■  vr  -  uq  +  a  +  (XG/m) 
v  ■  -ur  +  b  +  (YQ/m)  ) 

w  ■  uq  +  c  +  (ZG/m)  ^ 


(1-21) 
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(1)  The  primed  coordinate 
frame  is  fixed  in  the 
body . 

(2)  The  origins  are  co¬ 
incident  . 

(3)  The  Ox  and  Ox'  axes 
are  coincident. 


*  ■  P’  -  P 

♦  -  <K0)  +  ft  (p'-p)dt 

0 


a  -  a' 

o 

b  ■  b '  cos<f>  -  c '  s in4> 
o  o 

c  “  b^  sin4>  +  c^  cosb 

q  -  q'cosifi  -  r'sin<|! 
r  ■  q 1  sintj)  +  r '  costts 


the  Primed  and  Unprimed 


p 


0 


PB  -  L/A 

q  ■  (M  -  Ap’r)/B 
r  -  (N  +  Ap'q)/B 


(1-22) 


In  order 
aeroballistlc 


to  associate  the  unprimed  coordinate  frame  with  the 
frame,  first  the  spin  angular  momentum  is  defined  as 


^1  "  ABPB  +  AFPF  (1-23) 

It  is  assumed  that  the  fin-frame  will  have  axial  symmetry,  just  as  the 
body  does,  and  that  its  products  of  inertia  are  zero.  Equating  the 
rate-of-change  of  the  spin  momentum  to  a  total  rolling  moment, 


hi  -  AgPg  +  AyPy 


(1-24) 


However,  another  physical  variable  is  required  and  this  is  the  torque, 
T,  between  the  body  frame  and  the  fin-frame.  Let 


where 


A^P-c,  -  L,  +  T 


ABpB  "  LB  ”  T 


(1-25) 


Lj,  ■  aerodynamic  moment  applied  to  the  fin-frame 
Lg  •  aerodynamic  moment  applied  to  the  body-frame. 

It  can  be  seen  that  equations  (1-25)  satisfy  equation  (1-24) ,  provided  that 


L„  + 


4 


(1-26) 
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The  equations  which  have  been  developed  here  are  summarized  in 
Table  1.  Three  sets  of  equations  have  been  identified: 

(1)  Six  degrees-of-f reedom  expressed  in  body  axes 

(2)  Six  degrees-of-freedom  expressed  in  aeroballistic 

(pseudo-stability)  axes 

(3)  Seven  degrees-of-freedom  expressed  in  aeroballistic  axes 

The  first  set  is  most  general  but  it  presents  a  number  of  algebraic 
and  computational  difficulties.  The  second  set  introduces  the 
constraint  of  axial  symmetry  and  a  number  of  approximations  which 
are  convention  for  aeroballistics .  However,  it  does  not  include  an 
explicit  description  of  the  stabilized  and  controllable  fin-frame. 

This  second  set  has  been  included  here  in  order  to  show  a  relationship 
to  conventional  aeroballistic  theory.  The  third  set  of  equations 
introduces  a  seventh  degree-of -freedom:  the  spin  of  the  fin-frame 
independent  of  the  body.  This  seven  degree-of-f reedom  model  will  be 
the  basis  for  most  of  the  theoretical  design  considerations  for 
controlled  spinning  projectiles. 
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TABLE  1-1.  ALTERNATIVE  FORMULATIONS  OF  EQUATIONS  OF  MOTION 
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SECTION  2 


COORDINATE  TRANSFORMATIONS 


Three  dimensional  Cartiesian  vectors  transform  according  to  the 


(2-1) 


The  inverse  transformation  is: 


(2-2) 


and  the  transpose  of  the  matrix  is  also  its  Inverse.  The  matrix 
elements,  ,  are  called  "direction  cosines"  and  are  capable  of 
physical  interpretation. 


The  properties  of  these  coordinate  transformations  are  most 
familiar  to  engineers  in  terms  of  Euler  angles  (gimbal  angles). 

Such  a  presentation  is  given  in  Charters  [1] ,  and  none  will  be 
repeated  here.  A  less  familiar  treatment  of  the  coordinate  trans¬ 
formations  in  terms  of  the  components  of  a  rotation  vector  is  given 
in  Webster  [4]  and  will  be  discussed  here. 


The  rates-of-change  of  the  direction  cosines  depend  on  the 
present  values  of  the  direction  cosines  and  components  of  angular 
velocity.  The  complete  relationships  are: 


i 
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I 


11 

C21r 

c31q 

12 

m 

C22r 

- 

C32q 

13 

m 

c23r 

— 

c33q 

21 

m 

C31P 

- 

Cllr 

22 

m 

c32p 

- 

c12r 

23 

m 

c33p 

- 

c13r 

;31 

m 

cllq 

- 

°21p 

32 

m 

c12q 

- 

c22p 

l33 

- 

c13q 

- 

°23P 

(2-3) 


Each  rate-of-change  must  be  Integrated: 


5ij  "  ciJ(0)  +  J*  0 


ij 


dt 


(2-4) 


The  physical  significance  of  the  direction  cosines  can  be 
demonstrated  as  follows,  Consider  a  unit  vector  in  the  direction 
of  the  x°-axis.  It  can  be  expressed  as  a  column  vector: 

LSI 


l 

o 

o 


This  vector  will  transform  onto  the  unprimed  coordinate  frame  in 
accordance  with  equation  (2-1). 


(2-5) 
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Thus ,  C11  is  the  cosine  of  the  angle  between  the  x°-axis  and  the  x- 
axis;  c21  the  cosine  of  the  angle  between  the  x°-axis  and  the  y-axis; 
and  cai  is  the  cosine  of  the  angle  between  the  x°-axis  and  z-axis. 


The  argument  can  be  repeated  using  vectors 


0  "a  unit  vector  along  the  y°-axis 
1 
0 


■  a  unit  vector  along  the  z°-axis 


in  order  to  interpret  the  significance  of  the  remaining  . 

We  have  chosen  a  mnemonic  notation  for  the  vector  components 
because  we  have  found  that  it  facilitates  learning  and  remembering 
the  notation.  However,  there  is  some  merit  in  a  pure  index  notation 
and  such  a  notation  will  facilitate  the  present  discussion. 

The  equivalence  of  a  mnemonic  notation  and  an  index  notation  for 
a  right-handed  coordinate  frame  is  shown  in  Figure  (2-1).  The  axis 
are  assocated  with  the  indices  1,  2,  and  3  and  the  same  symbol  is 
employed  for  each  axis  or  for  each  vector  component.  The  same 
symbols  are  employed  for  each  unit  vector  and  the  unit  vectors  are 
distinguished  from  one  another  by  the  numerical  indices. 
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R  •  ix  +  jy  +  kz 

H  -  ip  +  jq  +  kr 
(ijk  are  unit  vectors) 


“1 


X  -  xl6l  +  x2e2  +  x3c3 

n  -  WjCj.  +  w2e2  +  0)363 

^ele2c3  are  unit  vectors) 


MNEMONIC  NOTATION 


INDEX  NOTATION 


Figure  (2-1).  Equivalence  of  Mnemonic  Notation  and  Index  Notation 
for  a  Right-Handed  Coordinate  Frame 


The  rates-of-change  of  the  direction  cosines  can  be  expressed 
in  the  index  notation: 


11 

= 

c2l“3 

c31Ll)2 

12 

a 

C22“3 

- 

c32u2 

13 

a 

C33w3 

— 

c33w2 

'21 

a 

c31w2 

- 

cllul 

22 

a 

c32u2 

- 

C12ul 

23 

a 

c33w2 

— 

c13ul 

:31 

a 

Cllai2 

- 

c21wl 

'32 

m 

C12U)2 

- 

c22ull 

:33 

m 

C13a)2 

- 

C23ul 

(2-7) 


Thus,  the  general  expression  is: 


ck+l,iwk-l  "  ck-l, iuk+l 


(2-8) 


where  k+1  and  k-1  are  to  be  evaluated  with  respect  to  the  following 
truth  table 


The  physical  basis  for  equations  (2-3)  or  (2-7)  can  be  demonstrated 
by  considering  small  rotations  instead  of  rates-of-change .  That  is 
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(2-9) 


$  -  /  pdt  -  /oj1dt 
6  =  Jqdt  *  /a)2^t 
^  =  /  rdt  »  J  Wgdt 


The  effect  of  these  small  rotations  is  to  change  the  projections  of 
a  unit  vector  in  the  inertial  frame,  x°  -  1,  on  the  axes  of  the  un¬ 
primed  frame.  The  effect  of  these  changes  is  depicted  in  Figure  2-2. 


The  Unprimed  Frame  Rotates  Through  Small  Angles 


Figure  2-2.  Physical  Basis  for  the  Change  in  Direction 
Cosines  Caused  by  Small  Rotations. 
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and  their  magnitude  can  be  written  by  inspection.  Thus 


AClj  "  C3j 6 
AC2j  "  C3j4>  "  Clj^ 
AC3j  “  Clj0  "  C2j* 


(2-10) 


The  several  vectors  described  in  Equations  1-3  can  be  transformed 
from  one  coordinate  frame  to  another  according  to  the  following 
schemes : 


T  . 


all  a12  a13 


a21  a22  a23  v  ; 


a31  a32  233  w 


°  ;  tJT  -  A  TJ° 


(2-11) 


bll  b12  bl3  u 


b31  b32  b33  w 


b21  b22  b23  I  I  v°  ’  B 


(2-12) 


C11  c12  c13  u 

C21  c22  c23  v°  ;  IT  -  C  U° 


C31  c32  c33  w 


(2-13) 


U"  dU  d12  d13  U 

d21  d22  d23  v  “  D  0 


d31  d32  d33 


(2-14) 
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SECTION  3 


THE  SEVEN  DEGREE-OF-FREEDOM  MODEL 


The  seven  degree-of-freedom  model  for  a  controlled  spinning 
projectile  has  been  derived  in  Section  1.  However,  it  will  be  helpful 
to  restate  the  equations  here: 


a  ■  X/m 
b  ■  Y/m 
c  ■  Z/m 


Accelerometer 

Equations 


(3-1) 


u  «  vr  -  wq 
v  ■  -ur  +  b 
w  ■  uq  +  c 


+  a  +  <XG/m) 

+  <Vm) 

+  (ZG/m) 


Force 

Equal lbruim 
Equations 


(3-2) 


hl  “  ABpB  +  afpf 
h2  ■  Bq 

h3  -  Br 


Angular 

Momentum 

Equations 


(3-3) 


PB  “  ^LB  "  T)/AB 
Pp  "  (Lp  +  T)/Ap 

q  -  (M  -  hjO/B 
r  «  (N  +  hxq)/B 


Moment 

Equilibrium 

Equations 


(3-4) 


These  equations  are  expressed  relative  to  the  aeroballistic  frame. 
The  seven  degrees  of  freedom  are  u,  v,  w,  Pp,  p^,  q  and  r. 


A  number  of  details  about  the  seven  degree-of-freedom  model 
must  be  clarified.  These  include: 
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(1)  Transformation  of  the  velocity  components  to  the  inertial 
frame  from  the  aeroballistic  frame. 

(2)  Transformation  of  the  gravity  components  to  the  aeroballis¬ 
tic  frame. 

(3)  Evaluation  of  the  aerodynamic  forces  and  moments  (X,  Y,  Z, 
Lf,  Lb,  m,  N). 


Equations  (3-1)  through  (3-4)  describe  the  motion  of  a  controlled 
spinning  projectile  with  respect  to  the  aeroballistic  frame.  The 
orientation  of  the  •aeroballistic  f  rame  in  inertial  space  is  determined 
by  the  angular  velocity  vector  (H)  where 


n 


r 


(3-5) 


These  vector  components  are  used  to  generate  the  transformation 
matrix  as  follows: 


C11  "  c21r  "  =31q 
c12  "  c22r  "  =  32q 
C13  "  C23r  "  =  33q 


C21  "  _cllr 
=22  "  “c12r 
623  "  “=13r 


=31  "  CUq 
=  32  "  c12q 
=33  "  =13q 


(3-0) 
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and  these  derivatives  are  integrated  as  follows 


cu  '  cij( 


r  ^ 

o)  +y  i 


C1J  dt 


(3-7) 


The  transformation  matrix  is  used  to  transform  the  u,  v,  w 
components  onto  the  inertial  frame: 


TT°  -C"1  0 


C11  C21  C31 


V  “  c12  c22  c32 

W°  c13  c23  C33 


(3-8) 


Components  of  wind  aloft  may  be  transformed  onto  the  aeroballistic 
frame  as  follows: 


CV 


cii  ci2  ci3  Nr 


21  22  23  VW 

i 

31  C32  c33  Nr 


(3-0) 


Also,  the  weight  of  the  projectile  can  be  transformed  from  the  inertial 
frame  to  the  aeroballistic  frame 


C11  c12  c13  0 
C21  c22  c23  0 
°31  c32  c33  W 


(3-10) 


I 


It  will  be  convenient  to  divide  equation  (3-10)  by  m  and  to  let 


g  ■  W/m 
gx  -  XQ/m 

Sy  "  Vm 

g2  -  ZQ/m 


(3-11) 


then 

^  -s 

*x  cn  ci2  cia  o 

Ky  “  C21  c22  c23  0  (3”^2> 

SZ  C31  c32  c33  8 

mJ  W  —  >-  - 

Consistent  with  equation  (3-12)  it  will  be  convenient  to  rewrite 
equation  (3-2)  as: 


u  -  vr  -  wq  +  a  +  gx 

v  -  -ur  +  b  +  g„ 

y 

w  ■  uq  +  c  +  gz 


(3-13) 


The  velocity  components,  u°,  v°,  w°, 
position  coordinates  as  follows: 


x°  -  x°(0)  +  ftu°dt 

</ 

y°  -  y°(0 )  +  ftv°dt 
0J 

z° '  »  z°(0  )  +  f1'w°dt 
0J 


are  integrated  to  obtain 


(3-14) 
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The  aerodynamic  forces  and  moments  are  expressed  in  the  following 


forms 

u -£^C 

where 


F  -  force  component 

M  -  moment  component 

C  -  dimensionless  coefficient 

9 

ipV  ■  dynamic  pressure 
S  -  reference  area 
d  •  reference  length 


The  complete  sets  of  forces  and  moments  which  are  required  are 


X 

- 

(pV2S/2)Cx 

Y 

SB 

( PV2S/2 )Cy 

> 

Z 

m 

(pV2S/2)Cz 

> 

lf 

m 

(pV2Sd/2)CLF 

lb 

m 

(pV2Sd/2)CLB 

M 

m 

(pV2Sd/2)CM 

> 

N 

m 

(pV2Sd/2)CN 

J 

(3-15) 


(.  3-16) 


I 


where  X,  Y,  Z  are  forces  in  the  x,  y,  z  directions,  respectively,  and 
3  •  L,  SI,  N  are  moments  about  the  x,  y,  z  directions,  respectively. 

The  moments  are  described  by  the  adjectives  rolling,  pitching,  and 
yawing,  respectively,  but  no  similar  adjectives  have  been  developed 
|  for  the  forces. 


u 


The  subscripts,  (  )B  and  (  )p,  have  been  introduced  to  distin¬ 
guish  between  the  rolling  moments  which  acts  on  the  body  (B)  and  the 
rolling  moments  which  acts  on  the  fin  frame  (F).  The  corresponding 
moment  coefficients  are  C^B  and  CLp,  respectively. 


It  became  common  practice  to  expand  the  dimensionless  aerodynamic 
coefficients  as  a  series  expressed  in  terms  of  a  number  of  dimension¬ 
less  variables.  The  more  important  dimensionless  variables  are: 


N 


a 

0 

(PBd/2V) 

<PFd/2V) 

(qd/2V ) 
(rd/2V) 

6 

A 


angle-of-attack  (rad) 
side-slip  angle  (rad) 
dimensionless  roll-rate  of  the  body 
(rad/sec) 

dimensionless  roll-rate  of  the  fin  frame, 
(rad/ sec) 

dimensionless  pitch-rate 
dimensionless  yaw-rate 
control  deflection  (rad) 
cant  angle  (rad) 


The  coefficients  of  each  series  tire  partial  derivatives  of  the  dimen¬ 
sionless  coefficients  with  respect  to  the  dimensionless  variables. 


However,  it  has  not  become  common  practice  to  employ  a  partial 
derivative  notation.  Instead,  a  notation  scheme  with  subscripts  is 
employed.  The  following  conventions  will  be  adopted  for  specifying 
§.  the  non-dimensional  aerodynamic  coefficients  ([1],  p.23). 
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CY  "  CYg  S  +  CYS  5Z 


cz  - 

-C_  a  - 

Za 

CZ6  6Y 

:L F  " 

+  CL  AF  +  cIe6 

3  ^Y  *  CIa6  a  6Z 

'LB  ” 

+ 

CM  “ 

C  a  +  i 
ma 

Cmq  (fv)  +  CmSp 

e  (w)  +  Cm6  6Y 

CN  “ 

c«6  8+0mr  (#)  +  Cnap  « 

(tt)  +  Cn5  5Z 

\  (3-17) 


J 


The  coefficients  of  the  dimensionless  terms  on  the  right  hand  sides 
of  the  equations  are  called  aerodynamic  derivatives  and  are  dimension¬ 
less.  The  following  relationships  are  valid  because  of  rotational 
symmetry. 


CY0 

"CZa 

CY5 

m 

CZ5 

C  « 

m 

-C 

n$ 

ma 

c 

m 

C 

nr 

mq 

t 

'nap 

m 

Cm8p 

Cn6 

m 

Cm6 

,F 

m 

CF 

^lafi 

(3-18) 


Thus,  a  total  of  19  aerodynamic  derivatives  of  which  12  are  independent 
will  be  used.  The  independent  derivatives  are  listed  below,  along 
with  an  alternative  notation. 
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r 


■ 


E 


p. 


r 

► 


► 


c 

X 

CX 

Cz, 

CZA 

Czc 

CZA 

^ma 

CMA 

^mq 

CMQ 

Cnap 

CNAP 

Cm5 

CMA 

CF 

Lip 

CLPF 

CF 

CLAF 

CB 

Lip 

CLPB 

CF 

'“lies 

CLBF 

CB  A 

CLAB 

axial  force 

normal  force  due  to  angle-of -attack 

normal  force  due  to  control  deflection 

moment  due  to  angle-of-attack 

damping  moment 

magnus  pitching  moment 

moment  due  to  control  deflection 

fin  rolling  moment  due  to  rate-of-roll 

fin  rolling  moment  due  to  cant  angle 

body  rolling  moment  due  to  rate-of-roll 

dihedral  effect 

body  rolling  moment  due  to  cant  angle 


The  symbols  in  tho  second  column  are  mnemonics  which  may  be  useful  in 
labelling  computer  printouts.  (The  derivatives  have  been  listed  here 
in  the  order  in  which  they  are  introduced  in  simulation  modules.) 

The  expressions  for  the  aerodynamic  derivatives  become: 


C  -  -C 
x 


cy  ' 
c_  - 


'LF 


'A 

“Cza  3  +  Cz6  6z 
"Cza  01  "  Cz6  ^y 


"LB 

CM 


C?A  AB 


3  5y  +  CA06  a{: 


m6  3y 


CN  "  “C 


:mq  (fv)  +  Cnap  3  +  C 

:ma  3  +  Cmq  l|v)  +  Cnap  a  +  Cm6  6 


(3-19) 


The  dimensionless  aerodynamic  derivatives  are  functions  of  Mach 
Number,  M.  For  the  present  these  functions  will  be  treated  as  second 
degree  polynomials,  e.g., 


Cx  -  aQ  +  a1M  +  a2M* 


(3-20) 


The  following  equations,  express  V,  a  and  8  in  terms  of  the  state 
variables  u,  v,  and  w,  and  the  wind's  components  u  ,  v  and  w  . 

w  w  w 

v2  -  (U  +  uw)2  +  (V  +  vw)2  +  (w  +  Ww)2  1 


8  ■  tan 


-  v  +  v 
-1  _ _w 

u  +  uZ 


(3-21) 


w  +  w_ 


oi  -  tan 


<u  +  uw)z  +  (v  +  vw)4 


The  most  important  equations  for  the  seven  degree-of -freedom 
model  are  summarized  In  Table  3-1. 


Fin  deflection  signs  refer  to  those  panels  whose  rotation  axis 
when  aligned  with  the  positive  y  or  z  axis  and  are  consistent  with 
positive  rotations  about  these  axes.  This  applies  to  either  panels 
which  deflect  together  (6)  or  panels  which  deflect  in  opposition  (A) 
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TABLE  3-1 

SUMMARY  OF  EQUATIONS  FOR  THE  SEVEN  DEGREE-OF-FREEDOM  MODEL 


a 

X/m 

"I 

u 

Y/m 

1 

Accelerometer 

0 

Equations 

c 

Z/m 

J 

\ 

vr  -  wq  +  a  +  g„ 

Force 

V 

-ur  +  b  +  g 

Equilibrium 

y 

Equations 

w 

uq  +  c  +  g2 

> 

hl 

ABFB  +  AFpF 

s  Angular 

’  Momentum 
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* 
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J 
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(pV2S/2)Cx 
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(pV2S/2)Cy 
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(pV2S/2)Cz 
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lf 

CPV2Sd/2)CLF 

Force  &  Moment 
/  Equations 

lb 
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TABLE  3-1 

SUMMARY  OF  EQUATIONS  FOR  THE  SEVEN  DEGREE -OF -FREEDOM  MODEL  (Continued) 


-C 


-c  e  +  c  ,  <$ 

za  zo  z 

-C  a  -  C  «  6 
za  z<5  y 


'n 


'LF 


'LB 


'M 


'N 


Clp  \W)  +  CL  AF  +  CIsi  6  6y  +  CJL65  ”’z 


a<5. 


-  c? 


© 


‘  =!.  aB 


^Aerodynamic 
i  Coefficients 


C  a  +  C 
ma  mq 


(tv)  +  Cnap  0 


&*) 


+  Cm6  5y 


Cma  0  +  Cmq  ($J  +  CQop  a  (^-)  +  Cm6  5 


This  completes  the  definitions  of  the  details  of  a  seven  degree- 
of-freedom  model  of  a  controlled  spinning  projectile.  All  of  this 
detail  can  be  embodied  in  a  single  software  module  which  can  be 
represented  by  the  block  diagram  shown  in  Figure  3-1.  This  can  be 
done  by  partitioning  the  model  into  the  following  parts: 

(1)  Acceleration  due  to  gravity 

(2)  Non-dimensional  terms 

(3)  Atmosphere 

(4)  Aerodynamic  forces  And  moments 

(5)  Equations  of  motion 

(6)  Coordinate  transformation  and  intergration 

(7)  Direction  cosines 

This  realization  is  described  in  the  appendix  as  the  realization  of 
Module  Number  60000  "Controlled  Spinning  Projectile".  Tho  appendix 
also  describes  the  realization  of  Module  Number  69220  "7-DOF  Projectile 
Including  Seeker  and  Canard  Deflection". 
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CONTROLLED 

SPINNING 

P_ROJECTILE 


Figure  3- 


Block  Diagram  Representation  of  a  Seven 
Degree-of-Freedom  Model  of  a  Controlled 
Spinning  Projectile. 


SECTION  4 


THE  FIVE-DEGREE  OF  FREEDOM  MODEL 

There  are  at  least  two  ways  of  simplifying  a  mathematical  model. 
One  way  is  by  an  arbitrary  introduction  of  simplifying  approximations. 
Another  is  by  a  careful  application  of  perturbation  theory.  Both 
ways  will  be  developed  here. 

The  first  way  (arbitrary  introduction  of  simplifying  assumptions) 
can  be  approached  as  follows: 

(1)  Assume  time  histories  or  constant  magnitudes  for  u  .and  for 
Pg  and  ignore  the  state  equations  for  u  and  p^. 

(2)  Evaluate  the  remaining  equations  using  the  assumed  magnitudes 
for  u  and  pg. 

(3)  Discard  terms  whose  magnitudes  are  judged  small  enough 
to  be  ignorable. 

This  technique  can  be  pushed  so  as  to  obtain  linear  equations,  but 
it  is  not  necessary  to  carry  it  that  far. 

The  equations  for  the  five  degree-of-f reedom  model  can  be  obtained 
by  deleting  the  appropriate  equations  from  (3-1),  (3-2),  (3-3)  and 
(3-4)  which  define  the  seven  degree-of-f reedom  model.  Thus 

Accelerometer 

Equations  (4-1) 


Force 

Equilibrium  (4-2) 

Equations 


b  -  Y/m 
c  ■  Z/m 

v  ■  -ur  +  b  +  g 
w  ■  uq  +  c  +  g_ 
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hl  "  AbPb  +  AFpF 

h2  -  Bq  l 

h3  =  Br 

✓ 

\ 

pp  -  ( L+T ) / Ap 

q  -  (M-h^/B  > 

f  -  (N+hiq)/B 


Angular 

Momentum 

Equations 


Moment 

Equilibrium 

Equations 


(4-3) 


(4-4) 


Equations  (4-1)  through  (4-4)  are  not  linear.  However,  they  can 
be  made  to  approach  linearity  by  the  following  approximations: 


h3PF  -  BrPp  *  0  > 

hgPp  -  BqPp  *  0 

hlr  “  (ABPB  +  AFPF)r  -  (ABPB)r  f 

hjQ  *  (abpb  +  ^pPp)q  **  ^  abpb^  ^ 

Substituting  (4-5)  in  equations  (4-1)  through  (4-4)  gives: 


(4-5) 


b  “  Y/m 
c  ■  Z/m 

v  ■  -ur  +  b  +  g 

6y 

w  -  uq  +  o  +  gz 
Pp  -  ( L+T ) / Ap 

q  •  (M-ABp^r)/B 
r  -  (N+ABPBq)/B 


*\ 


(4-6) 


It  remains  to  treat  the  force  and  moment  terms  (Y,  Z,  L,  M,  N)  in  a 
consistent  manner, 
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Note  that  the  terms  rpF  and  qpF  have  been  dropped  because  they 
are  second  order  quantities.  The  variables  v,  w,  p^,  q,  and  r  are 
expected  to  be  small,  so  their  products  are  second  order  quantities. 
The  product  terms  have  not  been  dropped  because  pF  ■  0.  The  magnitude 
of  pF  must  be  evaluated  and  the  consequences  of  its  being  non-zero 
must  be  included. 


The  required  force  and  moment  equations  can  be  selected  from 
equations  (3-15),  (3-16)  and  (3-17). 
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N 
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m  _C  6  +  C  „  6 
za  zf.  z 


C  a  -  C„  .fi  . 
zot  z  6  y 


Cj  -  C 


+  cIa  aF  +  c466  s6y  +  CJIBS  a6z 


“M 


C  a  +  C 
ma  mq 


+  c  B 

nap 


-,F 

+  C 


n<5  fiy 


CN  -Cma  *>  +  Cmq  (w)  *  Cnap  0  (tRt)  +  5: 


(4-7) 


Equations  (4-7)  can  be  used  to  express  the  forces  and  moments  in  the 
following  form  (see  equations  3-21). 
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Y  -  H  (V  +  vw)  +  ||  6; 


Z  - 


M 


N 


H  <w  +  V  +  H  s: 

pf  +  ti;  af 


IS  <w  +  ww)  +  If  q  +  If  (V  +  vw)  +  If  6. 

3N 


3q 

3N 


If  <v  ♦  vw>  ♦  %  r  *  g  (.  ♦  V  *  H  «. 


(4-8) 


01  course,  these  partial  derivatives  are  Q$Jt  non-dimensional, 


Equations  (4-8)  have  been  obtained  by  neglecting  the  product  terms 
6<5y  and  a$z  and  by  employing  the  following  approximations  (see 
equations  (3-21). 
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v  +  vw 

- v - 


w  +  w 
a  a  - ^ — 


>  (4-9) 
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The  expressions  for  evaluating  the  partial  derivatives  are: 
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Substituting  (4-8)  into  (4-0)  gives 


b  -  [ ( 3Y/3v)(v  +  vw)  +  ( 3Y/35 )Sy ] /m 


c  -  [ ( 3Z/3w)(w  +  ww)  +  ( 3Z/ 36 ) 5  ] /m 


v  •  -ur  +  b  +  g 
w  ■  uq  +  c  +  gz 

PF  -  U3L/3pF)PF  +  (3L/3Af)Af  +  Tl/Ay 
q  -  [(3M/3w)(w+ww)+(3M/3q)q+(3M/3v)(v+vw) 
,  +  ( 3M/36 )6y  -  ABp^r]/B 

r  -  [(3N/3v)(v+vw)+(3N/3r)r+(3N/3w)(w+ww) 
+  (3N/3S)6Z  +  ABp^q]/B 


(4- ID) 


>  (4-11) 
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We  introduce  some  more  notation 
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)  (4-12) 


(4-13) 


ik 


The  equation  for  is  quite  independent  of  the  other  equations. 
It  will  be  convenient  to  write  our  equations  in  the  following  form 


Pp  PF  +  PA  AF  +  PT  T 


■  >  — V  f  — 

v  Vv  0  0  Vr  v  0  V6  1  0  Vv  0 

w  -  0  Ww  Wq  0  w  +  0  010  Ww 

y 

q  Qv  Q^.  Qq  Qr  q  Qfi  0  0  0  Qy 

r  Rv  Rw  Rq  Rrj  r  ®  R<s  0  ®  Rv  Rw 

J  L  -K  J  l 


r  R. 


0  v6  1  0  yv  0 


+  V'6  0  0  1  0  Ww 

y 


0  R6  0  0  Rv  Rw 


(4-14) 


(4-15) 


b  Vv  0  0  0 

C  0  WOO 

m  W 

a  0  1/V  0  0 

6  1/V  0  0  0 


V,  0  0  V 

6  V 

z 


W.  0 

+  5y 
0  0 


0  0  0 


000  1/V 


0  0  0  1/V  0 


<4-fc  , 


The  five  degree-of-freedom  model  can  be  represented  by  a  block 
diagram  such  as  that  shown  in  Figure  4-1. 


Equations  (4-14)  (4-15)  and  (4-16)  are  written  with  respect  to  the 
aeroballistic  coordinate  frame.  The  control  deflections,  6^,  and  <$2, 
must  be  generated  in  this  frame  or  transformed  from  the  fin  frame. 
However,  the  gravity  components  and  the  wind  components  must  be 
generated  in  inertial  space  and  transformed  onto  the  unprimed  frame. 
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5  DOF  MODEL  OF  A 

CONTROLLED 

PROJECTILE 


Figure  4-1.  Block  Diagram  Representation  of  a  Five 
Degree-of-Freedom  Model  of  a  Controlled 
Projectile . 
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The  details  of  the  five  degree-of-freedom  model  can  be  realized 
most  easily  by  realizing  a  four  degree-of-freedom  model  first.  The 
next  section  (section  5)  should  be  read  at  this  point  and  forms  the 
basis  for  the  discussion  which  follows. 

The  realization  of  a  five  degree-of-freedom  time-variant  model  of 
a  controlled  projectile  is  shown  in  Figure  4-2.  This  model  is  a 
combination  of  the  modules  which  comprise  the  four  degree-of-freedom 
model  plus  additional  modules  to  simulate  the  moment  equilibrium 
about  the  x-axis.  The  spare  element  in  the  time-variant  element 
generator  can  be  used  to  generate  the  magnitude  of  p^.  Also, 
advantage  can  be  taken  of  the  fact  that  ■  u. 

A  five  degree-of-freedom  time-invariant  model  can  be  realized 
as  shown  in  Figure  4-3. 


4-0 


5UI 


Realization 


Figure  4—3.  Realization  of  a  5— DOF  Tiae— Invariant  Model 


SECTION  5 


THE  FOUR  DEGREE-OF-FREEDOM  MODEL 


The  four  degree-of-f reedom  model  can  be  obtained  readily  from 
the  five  degree-of-freedom  model  by  setting  pp  •  p^  •  0.  Then 


b  -  Y/m 
c  ■  2/m 


Accelerometer 

Equations  (5-1) 


v  ■  -ur  +  b  +  a  1  _ 

■y  \  Force 

w  -  ua  +  c  +  e  I  Equilibri 

w  uq  +  c  +  g2  j  Equations 


(5-2) 


hi  "  Vb 

h2  ■  Bq 
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Angular 

Momentum 

Equations 


(5-3) 


q  ■  (M  -  h1r)/B 
f  ■  (N  +  hiq)/B 


Moment 

Equilibrium 


(5-4) 


These  equations  will  lead  to  equations  which  are  equivalent  to 
equations  (4-15),  (4-16)  and  a  subset  of  (4-13).  These  equations  are 
rewritten  here  so  as  to  have  a  compact  summary . 
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(5-6) 


(5-7) 


Equations  (6-6)  and  (5-7)  are  linear,  but  the  matrix  elements  may  be 
time  variant. 


In  order  to  accommodate  the  time-variant  nature  of  the  four 
degree-of-freedom  model  it  will  be  organized  as  shown  in_ Figures  5-1, 
5-2  and  5-3.  The 'time-varying  coefficients  are  the  magnitudes  given  in 
in  equation  4-5.  These  coefficients  can  be  determined  in  a  7  DOF  sim¬ 
ulation  and  then  expressed  as  polynominals .  Quadratic  polynomials 
will  be  reasonable  at  this  time. 
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Figure  5-1.  Block  Diagram  of  4-DOF  Force  Equilibrium  Equations. 


Figure  5-2.  4 -DOF  Moment  Equilibrium  Equations 


Figure  5-3.  Block  Diagram  of  a  4— DOF  Hodel  of  a  Controlled  Projectile 


Quadratic  polynomials  can  be  generated  by  module  number  8130. 

Ten  of  tbese  can  be  combined  in  a  single  module  as  is  shown  in  Figure 
5-4  . 


The  realizations  of  time-variant  and  time- invariant  4  degree-of- 
freedom  models  are  described  and  compared  in  Figure  5-5. 


Figure  5-4 .  Block  Diagram  of  Time  Variant  Elements  for  4 -DOF  Model . 


4  DOF  LINEAR  4  DOF  LINEAR 

MODEL  OF  A  MODEL  OF  A 


Comparison  of  a  Tine-Variant  and  a  Time- Invariant  4-DOF  Linear 
Model  of  a  Controlled  Projectile. 


SECTION  6 


SYSTEM  DESCRIPTIONS  EMPLOYING  EXACT  EXPRESSIONS  FOR  THE  COORDINATE 

TRANSFORMATIONS 

The  approximate  five  and  four  degree-of-freedom  models  can  be 
employed  in  conjunction  with  exact  or  approximate  expressions  for  the 
coordinate  transformations.  The  approximate  expressions  are  more 
eonomical  of  computing  power  and  permit  linearization  and,  therefore 
they  will  be  used  whenever  possible.  However,  the  exact  expressions 
will  be  useful  for  error  control. 

The  complete  description  of  the  projectile,  which  parallels -that 
presented  for  the  seven  degree-of-freedom  model  (see  Section  3)  is 
represented  in  Figure  0-1. 

The  target  can  be  represented  by  the  block  diagram  shown  in 
Figure  6-2. 

The  relative  motion  of  the  projectile  and  the  target  is  described 
by  the  following  equations: 
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8210 
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Then , 


RR  -  Ax°  Au°  +  Ay°  Av°  +  Az°  Aw° 
R  -  ^Au°  +  Avo  +  Awo 


(6-0) 


A  seeker  mounted  in  the  aeroballistic  frame  (unprimed  coordinate  frame) 
will  develop  boresight-error  angles  which  depend  on  the  position  of  the 
target  relative  to  the  projectile,  resolved  onto  the  unprimed 
coordinate  frame.  Thus, 


S3 T  -  C  AX°  (see  Appendix  A) 


The  boresight-error  angles  will  be 


(0-7) 


(0-8) 


The  boresight-error  angles  measure  the  relative  position  of  the 
target  and  projectile,  but  they  also  include  the  effect  of  attitude 
change  of  the  projectile.  This  is  not  obvious  from  the  exact 
equations;  but  it  will  be  obvious  from  the  approximate  equations  which 
are  developed  in  the  following  section. 


The  relative  motion  of  the  target  and  the  projectile  can  be 
represented  by  the  block  shown  in  Figure  6-3.  This  block  also 
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RELATIVE  MOTION 
AND 

STRAPPED  -  DOWN 
SEEKER 


Figure  8-3.  Block  Representing  the  Relative  Motion  of  Target  and 
Projectile  and  the  geometric  processing  performed  by 
a  Strapped-Down  Seeker. 
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» 


represents  the  geometric  processing  done  by  a  strapped  down  seeker 

and  generates  the  boresight-error  angles,  X  and  X  .  These  functions 

y  z 

are  valid  for  the  7-DOF,  5-DOF,  and  4-DOF  models. 
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SECTION  7  .  | 

SYSTEM  DESCRIPTIONS  EMPLOYING  APPROXIMATE  EXPRESSIONS  FOR 
THE  COORDINATE  TRANSFORMATIONS 

System  descriptions  employing  approximate  expressions  for1  the  > 

coordinate  transformations  permit  linearization  in  a  manner  comparable 
to  that  which  has  been  introduced  for  the  four  and  five  degree  of 
freedom  models.  Such  treatment  permits  the  linearization  of  the 
complete  interception  problem.  The  linearized  state-space  equations  > 

for  the  physical  plant  are: 

» 

(7-1) 

•  i 


> 


* 
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> 


*(t)  -  A(t)x(t)  +  B(t)u(t) 

£(t)  -  C(t )x(t)  +  D(t )u(t) 

where 

x(t)  -  a  state  vector 
£(t)  ■  an  output  vector 

u(t)  "  an  input  vector 

The  state  vector  will  include  the  following  terms: 

PF  angular  velocity  of  the  fin  frame  (rad/'sec) 

translational  velocity  (ft/sec) 

q  ) 

>  angular  velocity  (rad/sec) 


7-1 


6  rotations  (rad) 


translational  perturbations  (ft) 


<5  control  deflection  (rad) 


The  output  vector  will  include  the  following  terms: 


translational  acceleration  (ft/sec^) 


a  angle-of-attack  (rad) 
0  side-slip  angle  (rad) 


relative  position  components  (ft) 


y  )  look  angle  components  (rad) 


The  input  vector  will  include  the  following  terms: 


wind  components  (ft/sec)  Ap  fin  cant  angle  (rad) 


>  target  maneuvers  (rad/ sec) 
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The  approximate  expressions  developed  in  this  section  will 

employ  a  coordinate  frame  which  is  fixed  to  the  cruciform  canards 
F  F  F 

(the  x  ,  y  ,  z  ).  This  convention  makes  it  possible  to  include  the 
effects  of  small  perturbations,  $ ,  in  the  roll  angle  of  the  canard 
frame,  The  perturbation  equations  for  the  direction  cosines  are 
similar  to  those  which  have  been  developed  in  equations  (2-1Q). 
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The  resolution  of  gravity  onto  the  fin  frame  is  described  by: 
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The  usual  practice  will  be  to  include  the  effect  of  gravity  in  the 
reference  trajectory.  Then  it  is  only  necessary  to  treat  Ag  in  the 
linearized  equations.  Thus: 
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The  transformation  of  the  wind  components  onto  the  fin  frame  is 
given  by: 
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No  perturbation  terms  are  included  in  the  transformation 
because  the  wind  components  are  considered  to  be  small  quantities . 

The  velocity  components,  computed  in  the  fin  frame,  must  next 
be  resolved  onto  the  inertial  frame.  Thus: 
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Expanding  equation  (7-7)  and  remembering  that  Vp  and  w£,  are  small 
perturbations : 


(7-8) 
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The  target's  velocity  can  be  specified  so  that  -  u.  Then  the 

target's  velocity  can  be  resolved  onto  the  inertial  frame  in  a  similar 
manner  and: 
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The  position  coordinates  of  the  target  are: 
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and  the  coordinates  of  the  projectile  are: 
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The  relative  position  components  of  target  and  projectile,  defined 


inertial 

frame  are : 
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It  can  be  seen  that  the  components  of  A*X° ,  AY0,  A*Z°  are  functions 

of  time,  whereas  the  components  of  Ax° ,  A$°,  Az°  are  functions  of  the 
F  F 

perturbations  v£,  wp,  9  and  These  relative  displacement  compo¬ 
nents  must  now  be  transformed  back  to  the  fin  frame. 
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The  right  hand  side  of  equation  (7-15)  expands  into  four  terms  which 
are  treated  below,  one  at  a  time. 


The  first  term  is: 
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The  second  terra  is: 
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The  third  term  is: 
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The  fourth  term  is: 
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The  expansion  of  this  term  contains  only  products  of  perturbations, 
all  of  which  can  be  ignored.  Then  the  expansion  of  equation  (7-16) 
is: 
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Equation  (7-20)  describes  the  relative  position  of  the  target 
and  projectile  resolved  onto  the  fin-frame.  It  can  be  seen  that 
this  equation  encapsulates  a  great  deal  of  detail  which  would 
require  a  large  amount  of  computation  if  the  perturbation  equations 
had  not  been  developed. 
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It  may  be  useful  to  identify  the  terms  which  appear  on  the 
right-hand  side  of  equation  (7-20),  thus: 


~F  ~F  ~F 
nx  ,  Ay  ,  Lz 


are  relative  position  components  deter¬ 
mined  during  the  reference  trajectory. 


F 

«P 


projectiles  component  of  velocity  along 
the  x-axls  of  the  fin-frame. 


vp,  wp,  6,  ty 


perturbations  of  the  projectiles  motion 
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T  **p  v 

Thus,  Ax  ,  Ay  ,  Az  ,  Up,  are  time-variant  coefficients  to  he  deter¬ 
mined  from  the  computation  of  a  reference  trajectory  while  v^  wp, 

6 ,  w  are  state-variable  perturbations  to  be  determined  from  u  linear 
model  of  an  interception. 

The  equations  developed  in  this  section  are  an  alternative 
to  the  non-linear  equations  developed  in  sections  (1)  through  (3) 
and  in  Section  (6).  All  of  these  sections  were  developed  with  the 
description  of  motion  referred  to  the  aeroballistlc  coordinate  frame. 
This  convention  has  merit  in  facilitating  understanding  of  complex 
motions.  However,  the  aeroballistlc  frame  convention  does  not  permit 
a  roll-perturbation,  to  be  derived  as  a  small  quantity  of  the 
first  order.  Formulating  the  description  of  the  motion  with  respect 
to  the  fin  frame  permits  retaining  the  roll  perturbation  as  a  small 
quantity  of  the  first  order  in  a  linear  model. 


Equation  (7-20)  has  been  helpful  in  understanding  the  roll- 
control  of  the  fin-frame.  It  will  also  be  of  considerable  importance 
when  optimal  design  formulations  are  undertaken. 
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APPENDIX  A 


7-DOF  PROJECTILE  MODULE  DESIGN  DETAILS 


A . 1  INTRODUCTION 


This  appendix  describes  the  details  of  a  moculd  which  simulates 
seven  degrees-of-f reedom  of  a  CHAMP  projectile.  The  following 
modules  are  involved  in  this  production: 

8210  Acceleration  due  to  gravity 

57130  Non-dimensional  terms  with  wind 
0200  Atmosphere 

57110  Aerodynamic  derivatives 

60050  Aerodynamic  forces  and  moments 
57030  Equations  of  motion 

860  Direction  cosines 

870  Transformation  and  integration 
60060  Controlled  spinning  projectile 

60220  7-DOF  projectile  including  seeker  and  canard  deflection 

Details  of  these  modules  are  presented  in  the  following  sections  of 
this  appendix.  Background  material  on  the  modular  software  system 
has  been  published  in  the  CHAMP  Phase  I  Final  Report,  Appendix  B. 

A. 2  ACCELERATION  DUE  TO  GRAVITY  (Module  #8210) 


The  force  due  to  gravity  can  bo  expressed  conveniently  in  the 
inertial  coordinate  frame.  The  components  in  this  frame  are: 
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(where  g  ■  acceleration  due  to 
gravity) 
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These  components  can  be  transformed  to  the  unprimed  coordinate 
frame  by  the  following  transformation: 
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and  equation  (A-4)  gives  the  terms  which  are  required  in  the  force 
equilibrium  equations,  (3-2). 


This  computation  can  be  performed  by  Module  No.  8130  which  was 
coded  directly  in' MARCO-10.  The  module  is  represented  below  and 
is  summarized  in  Table  A-l. 


FORCE 
DUE  TO 
GRAVITY 


Figure  A-l. 


Module  No.  8210  Force  Due  to  Gravity. 
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TABLE  A-l 

MODULE  NO.  8210  -  ACCELERATION  DUE  TO  GRAVITY 


Algebraic 


L  M 

Symbol 

Dimensions 

Explanation 

Connections 

3  1 

c13 

dimensionless 

direction  cosine 

2 

c23 

dimensionless 

direction  cosine 

3 

c33 

dimensionless 

direction  cosine 

Parameters 

4  1 

mg 

lbs 

weight  of  the  projectile 

Intermediates 

7  1 

XG 

lbs 

force  components  due  to 

2 

yg 

lbs 

gravity 

3 

ZG 

lbs 

A .  3  NON-DIMENSIONAL  TERMS  (Modules  #57000  8t  #57130) 

The  Inputs  to  thin  module  are: 

z  -  the  negative  of  the  altitude,  h  (feet) 

u,v,w  ■  velocity  components  ( feet/second) 

pF,PB,q,r  ■  angular  velocity  components  (radians/second) 

The  functions  of  this  module  are: 

(1)  to  transform  z  into  h*-z,  which  is  required  as  an  input 
by  the  ATMOSPHERE  module  (#9200)  . 
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(2)  to  derive  V,  a,  and  B  from  u,v,w 

(3)  to  derive  the  non-dimensional  terms  Ppd/2v,  Pgd/2v, 
qd/2v,  and  rd/2v,  where  d  is  a  reference  length. 


The  composite  module  can  be  implemented  with  standard  modules  as 
shown  in  Figure  A-2. 


This  same  module  can  be  employed  to  evaluate  the  non-dimensional 
terms  required  to  describe  the  effect  of  winds,  provided  the 
following  replacements  are  made: 


u  «■  u+u^ 
v  *■  v+v 

w 

w  +■  w+w„ 
w 

where  u,v,w  are  state  variables  resolved  on  the  unprimed  coordinate 
frame.  uw>vw*ww  are  wiQd  components  resolved  on  the  unprimed 
coordinate  fra-T. 


Now  the  wind  components  are  best  described  in  the  inertial  frame, 
so  that : 
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The  transformation  and  additions  have  been  combined  in  a  single 
module  as  shown  in  Figure  A-3  and  Table  A-2. 
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Figure  A-3.  Module  No.  57130  Non-Dimensional  Terms  (with  wind) 


TABLE  A- 2 

MODULE  NO.  57130  NON-DIMENSIONAL  TERMS  (WITH  WIND)  DICTIONARY 


L  M  Algebraic 

-  -  Symbol  Dimensions  Explanations 

Connections 


3  1 

z 

ft 

z  position  coordinate 

2 

u  ] 

velocity  components  re¬ 

3 

V  \ 

ft/sec 

solved  onto  the  unprimed 

4 

W  | 

coordinate  frame 

5-13 

ciJ 

dimensionless 

direction  cosines 

14 

velocity  components  of  the 

15 

16 

v°  1 

vw  1 

w° 

ww 

ft/sec 

wind  resolved  onto  the 

inertial  coordinate  frame 

17 

Pp 

rad/ sec 

x-component  of  angular 
velocity  of  the  fin  frame 

18 

A 

PB 

rad/ sec 

x-component  of  angular 
velocity  of  the  body 

19 

20 

*  } 
r  1 

[  rad/sec 

1 

angular  velocity  components 

Parameters 

4  1 

d/2 

ft 

reference  length 

Intermediates 

7  1 

uw 

velocity  components  of  the 

2 

vw 

ft/sec 

wind  resolved  into  the 

unprimed  coordinate 

3 

ww 

frame 
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TABLE  A-2  (continued) 

MODULE  NO.  57130  NON-DIMENSIONAL  TERMS  (WITH  WIND)  DICTIONARY 


L  M  Algebraic 

—  —  Symbol  Dimensions  Explanation 

Intermediates 


4 

u+u 

w 

velocity  components  of 

5 

6 

v+vw 
w+ww  . 

« 

ft/sec 

the  relative  wind  resolved 

onto  the  unprimed  coord¬ 
inate  frame 

7 

h 

ft 

altitude 

9 

V 

ft/sec 

magnitude  of  the  relative 

wind 

11 

a 

rad 

angle-of-attack 

12 

6 

rad 

side  slip  angle 

13 

Ppd/2v  > 

dimensionless 

non-dimensional  terms 

14 

Pgd/2v 

dimensionless 

non-dimensional  terms 

15 

qd/2v 

f 

dimensionless 

non-dimensional  terms 

16 

rd/2v 

J 

dimensionless 

non-dimensional  terms 

> 


•  > 


> 
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A. 4  ATMOSPHERE  (MODULE  NO,  9200) 

The  relationships  necessary  to  define  atmospheric  properties 
as  a  function  of  altitude  are  based  on  standard  equations  of  fluid 
statics.  The  numbers  used  in  the  equations  are  those  of  the  1959 
ICAO  Standard  Atmosphere  and  its  extension.  The  Model  is  piecewise 
continuous  in  three  segments.  It  is  based  on  a  constant  lapse  rate 
from  sea  level  to  the  tropopause  (36,089  ft.).  The  atmosphere  is 
assumed  to  be  Isothermal  from  36,089  ft.  to  82,021  ft.  From  82,021 
ft.  to  154,199,  the  static  temperature  increases  at  a  constant  rate. 

The  outputs  of  the  atmosphere  model  are  uniquely  determined 
given  the  altitude  as  an  input.  These  outputs  include  static 
pressure  ratio,  and  local  speed  of  sound.  In  this  particular  formu¬ 
lation,  one  subsidiary  calculations  is  also  performed.  The  air¬ 
speed  is  used  as  an  input  to  the  model  so  that  Mach  No.  can  be 
calculated.  Thus,  having  specified  altitude  and  airspeed,  the  model 
determines  static  pressure  ratio  and  Mach  No.  for  use  in  calculating 
aerodynamic  forces  and  torques.  This  module  is  coded  in  MARCO-10. 

Definitions : 


a  is  local  speed  of  sound  in  ft/sec 

X  is  ratio  of  local  static  pressure  to  sea  level 

static  pressure 

M  is  Mach  Number 

tr  is  the  ratio  of  local  static  temperature  to  sea 

level  static  temperature 

V  is  true  airspeed  in  ft/sec 

h  is  altitude  above  sea  level  in  ft. 
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Connections : 


V 


h 


Parameters ;  None. 
Initial  States:  None. 


Intermediates : 

tr  -  l-( 6 . 8754  x  10‘6)h 
X  -  tr8’2861 


0  <  h 


tr  -  0.75187 


X 

tr 

X 


-5, 


1 . 20586  <~‘-8063  11  10  >h 


0.40100  +  3.1732  x  10~6h 


0.02450(0.85383  +  (4.2204 
x  lO-^h)-11*388 


36089 


82021 


tr  -  1.0 

.  (-.322004 +  (-4.115414  x  10_5)h) 

A  "*  © 


a  ■  1118 . 9  \/  tr 
M  -  V/a 


36080 


h  <  82021 


h  <  154190 


|  h  >  154190 


States:  None. 


A . 5  AERODYNAMIC  FORCES  AND  MOMENTS 


The  inputs  to  this  module  are: 

a.  The  control  deflections  (<5y  and  6g) 

b.  The  cant  angle  (Ap) 

c.  The  pressure  ratio  and  Mach  Number  (X  and  M) 

d.  The  non-dimensional  terms  (a,  B,  prd/2v,  p„d/2v,  qd/2v, 

rd/2v)  ?  ® 

The  most  important  outputs  are: 

a.  The  components  of  aerodynamic  force  (X,  Y,  Z) 

b.  The  components  of  aerodynamic  moment  (Lp,  LQ,  M,  N) 

The  simulation  is  synthesized  in  two  steps.  First,  the 
aerodynamic  derivatives  are  evaluated  as  functions  of  Mach  Number. 
This  is  formulated  in  Modules  Number  57110  and  8130.  Then  the 
derivatives  are  combined  to  evaluate  the  non-dimensional  force  and 
moment  coefficients  and  are  multiplied  by  the  appropriate  dimensional 
terms  to  yield  force  components  in  pounds  and  moment  components  in 
foot-pounds.  This  is  accomplished  in  Module  Number  69050. 

A. 5.1  AERODYNAMIC  DERIVATIVES  (Module  No.  57110) 

The  several  aerodynamic  derivatives  are  expressed  as  quadratic 
functions  of  Mach  Number;  i.e., 

C  -  AQ  +  A^M  +  AjM2 

where  the  coefficients  AqI  A^,  Ag  are  fitted  to  the  best  available 
wind  tunnel  data  or  estimates  in  an  independent  procedure.  The 
derivatives  which *are  employed  in  the  simulation  are: 
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axial  force 

Cr-  normal  force  due  to  angle  of  attack 

CZ(5  normal  force  due  to  fin  deflection 

C  pitching  moment  due  to  angle  of  attack 

ma 

C  damping  in  pitch 

mq 

C _  Magnus  moment 

nap 

Cm6  pitching  moment  due  to  fin  deflection 
F 

C'p  roll  damping  of  the  fin 

F 

C"^  rolling  moment  due  to  fin  cant 

roll  damping  of  the  body 
F 

dihedral  torque  on  the  canard  frame  due  to  fin 
deflection 

body  rolling  moment  due  to  cant  angle 

The  computation  also  evaluates  the  dynamic  pressure: 
q  -  1482.52M2X  (lbs/ ft2) 

These  terms  are  evaluated  in  Module  57110  presented  in  Figure  A-5. 

A. 5.2  AERODYNAMIC  FORCES  MOMENTS  (Module  No.  09050) 

Module  Number  09050  evaluates  the  dimensional  forces  and 
moments  which  act' on  the  projectile.  An  additional  aerodynamic 
term  is  required  beyond  those  included  in  Module  Number  57110  and 
this  is  added  using  Module  Number  8130.  Module  80050  is  illustrated 
in  F igure  A-0 . 
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Figure  4-5 .  Module  Muuber  57110  Aerodynamic  Derivatives 


A . 6  CONTROLLED  SPINNING  PROJECTILE  (Module  66060) 


Module  69060  simulates  the  seven  degrees-of-freeriom  of  a 
controlled  spinning  projectile  in  a  manner  which  is  consistent  with 
the  theory  derived  in  Section  3  of  the  main  body  of  thie  volume. 

The  coordinate  frame  for  the  module  ie  made  to  correspond  to 
the  pseudo-stability  axes  by  setting  p  ■  [09060,3,5]  equal  to  zero. 
The  module  accepts  wind  components  relative  to  the  inertial  coordin¬ 
ate  frame  and  it  computes  the  projectiles  motion  relative  to  the 
inertial  coordinate  frame. 

Attention  should  be  called  to  the  sign  conventions  for  wind. 

The  wind  components  have  been  taken  as  positive  when  they  are  in 
the  direction  of  the  negative  coordinate  axes.  This  permits  the 
total  relative  wind  components  to  be  written  as  u+uw,  v+vw,  w+ww  . 
consistent  with  equations  3-21. 
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A.  7 


7- DOF  PROJECTILE  INCLUDING  SEEKER  AND  CANARD  FRAME  (Module  69220) 


Module  69220  embodies  the  seven  degree-of-freedom  simulation  of 
Module  69060.  In  addition,  it  includes  the  following  features: 

(1)  It  computes  the  relative  position  of  the  target  and 
projectile  in  the  inertial  coordinate  frame. 

(2)  The  relative  position  is  resolved  onto  the  aero- 
ba.llistic  coordinate  frame  and  look-angle  components 
are  evaluated. 

(3)  The  angle,  between  the  fin-frame  and  the  body 
is  evaluated. 

(4)  The  angle,  is  used  to  resolve  the  canard  deflec- 

e  F 

tion  in  the  canard  frame,  SY>  onto  the  aeroballistic 
frame  and  the  look-angle  components  onto  the  fin- 
frame  . 


